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Abstract 

We investigate the nucleoli's transverse spin-dependent structure function 
g2{x,Q'^) in the framework of the operator product expansion and the renor- 
mahzation group. We construct the complete set of the twist-3 operators for 
the flavor singlet channel, and give the relations among them. We develop an 
efficient, co variant approach to derive the anomalous dimension matrix of the 
twist-3 singlet operators by computing the off-shell Green's functions. As an 
application, we investigate the renormalization mixing for the lowest moment 
case, including the operators proportional to the equations of motion as well 
as the "alien" operators which are not gauge-invariant. 
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Recently, the nucleoli's transverse spin-dependent structure function g2{x,Q'^) has 
been observed by measuring asymmetry in the deep inelastic scattering using the trans- 
versely polarized target. The structure function g2{x, Q"^) plays a distinguished role in spin 
physics because it contains information inaccessible by the more familiar spin structure 
function Q^) 0, 0]: g2{x,Q'^) is related to the nucleon's transverse polarization and 
to the twist-3 operators describing the quark-gluon and three gluon correlations in the 
nucleon. 

In the framework of the operator product expansion, not only the twist-2 operators but 
also the twist-3 operators contribute to g2{x,Q'^) in the leading order of l/Q^P, The 
Q^-evolution of g2{x, Q^) is governed by the anomalous dimensions, which are determined 
by the renormalization of these operators. A characteristic feature of the higher twist 
operators is the occurrence of the complicated operator mixing under the renormalization: 
Many gauge-invariant operators, the number of which increases with spin (moment of 
the structure function), mix with each other^j. Furthermore, the operators which are 
proportional to the equations of motion ("EOM operators"), as well as the ones which 
are gauge- variant ("alien operators"), also participate in the mixing]^]. 

There have been a lot of works on the Q^-evolution of g2{x,Q'^)- Most of them discussed 
the flavor non-singlet case^, ^j. Only a few works treated the singlet case0: Bukhvostov, 
Kuraev and Lipatov derived evolution equations for the twist-3 quasi-partonic operators. 
Recently, Miiller computed evolution kernel based on the nonlocal light-ray operator tech- 
nique, and obtained the identical results. However, both of these two works employ a 
similar framework based on the renormalization of the nonlocal operators in the (light- 
like) axial gauge. Balitsky and Braun also treated the nonlocal operators although they 
employed the background field method. On the other hand, a covariant approach based 

on the local composite operators is missing. Furthermore, some subtle infrared problem 

*Recently, it has been proved that the twist-3 non-singlet structure functions obey simple DGLAP 
evolution equation for Nc oopl|. 
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occurring in the renormalization of the generic flavor singlet operators has been empha- 
sized in Ref.0. Therefore, the computation of the anomalous dimensions in a covariant 
gauge in a fully consistent scheme is desirable and should provide a useful framework. 

We develop a covariant framework to investigate the flavor singlet part of g2{x,Q'^) 
based on the operator product expansion and the renormalization group, by extending 
our recent work on the flavor non-singlet part]^, ^. The purpose of this letter is twofold. 
Firstly, we provide all necessary stuffs and convenient techniques: We list all relevant 
twist-3 flavor singlet operators appearing in QCD. We discuss the relations satisfied by 
these operators]^. In particular, we obtain a new operator identity, which relates the 
gluon bilinear operator with the trilinear ones. Based on these developments, we give a 
basis of independent operators for the renormalization. We describe a general framework 
to perform the renormalization of the twist-3 flavor singlet operators in a covariant gauge 
by computing the off-shell Green's functions. An advantage of our approach is that an 
infrared cut-off is provided by the off-shellness of the external momenta. This allows us 
to assess unambiguously the renormalization constant ; otherwise the calculation might 
be subtle and potentially dangerous|Q. A convenient projection technique^, ||] is also 
introduced to simplify the calculation. Secondly, we apply our framework to the lowest 
(n = 3) moment case, and demonstrate the consistency and the efficiency of our method: 
We investigate in detail the mixing structure of the higher twist singlet operators under the 
renormalization. The EOM operators as well as the alien operators should be included as 
independent operators. We also clarify the connection of our approach with other works 
based on the on-shell calculation without the EOM operators nor the alien operators. As a 
byproduct, we confirm the prediction of the Q^-evolution obtained in previous works[^, |^, 
which will be subject to future experimental studies. 

We follow the convention of Refs.^, ^. We list the twist-3 flavor singlet operators 
which contribute to the moment Jq dxx^'^~^g2{x, Q"^) [n = 3, 5, 7, ■ ■ ■). Firstly, there is a 



2 



set of operators which are bihnear in the quark fields. These operators can be obtained 
from the non-singlet operators, eqs.(l)-(4) of Ref.0 by removing the fiavor matrices Aj. 
An independent basis can be chosen as (see Refs. ||, [K 



RZ"'-' = ^{Vi-Vn-i^i + Ui + Un-i-i) (/ = l,---,n-2) , (1) 
= i'^-^Smi^j^YD''' ■ ■ ■ D^^-^-i^^-^i) - traces , (2) 



n — 1 



2n 

+ ■ ■ ■ /^^'"-^^'^n-i [ijp- - traces , (3) 



where 



Vi = i'^gStp-f^Df^'-'-G"'" ■■■Df'^-^Y"~''4^ -traces, 
Ui = €'~^gSi)D^^ ■ ■ ■ G"'" ■ ■ ■ Df'^-^Y^-'tP - traces. 

In the above equations, (- traces) stands for the subtraction of the trace terms to make 
the operators traceless and is the covariant derivative. S means symmetrization over 
fii and g the QCD coupling constant, m represents the singlet component of the quark 
mass matrix. The operators in eq.(|l|) contain the gluon field strength G^,^ and the dual 
tensor G^u = \e^uai3G°'^ explicitly. 

Secondly, we have the following "new" operators (see e.g. Ref . Jim ) : 



traces, (4) 



_ ■n-2g ^ \^G''f^iD^2 . . . ^-/., . . . D^'-^Gf-'l - traces (5) 

n — 1, 

{1 = 2 



2 " 

_ ■ ■ ■ D^'^^-^Y{--d>'--'{d''Al) + ^r'"(5^"-^x'')x'} (6) 

a 

—traces, 

+ gi^t^j^-'ij + i^^-i {d'^AD - gr''%d^'"''x'')xl - traces. (7) 
a 



Here A antisymmetrizes cr and fii. The gluon field and the covariant derivative D'^ are 
in the adjoint representation, x is the ghost field, is the color matrix as [t", t*] = if^^H"^, 
Tr (rt^) = and a is the gauge parameter. T„ ^ is trilinear in the gluon field strength 
G^i, and the dual tensor, and thus represents the effect of the three gluon correlations. It 
satisfies the symmetry relation = T^'^_i^'^~^ . Tn,E is the gluon EOM operator; 

it vanishes by the use of the equations of motion for the gluon, although it is in general a 
nonzero operator due to quantum effects. Tn.B is the BRST-exact alien operator 0] which 
is the BRST variation of the operator: i'^-^S{G''^'^ D^'^ ■ ■ ■ D''"-^Yd^'"-'x^ - traces. 

We note that the gluon bilinear operator eq.(^) is related to the trilinear ones eq.(^ 

by 



n-l 



- {l-2)Cr' (n-/-2)Crf 



1=2 



2{n-l) 2(n-r 

yn-2 



2(n-2)(n-l)^ ' 

n~2 



1 \lrp(Jfli-- 

-'-nJ 



where = n\/[r\{n — r)!]. To derive this relation, we have used [D^, Dy] = —igG^^ and 

■D(jGi/a ~\~ -D pG c^(j -\- ID aG 0, D^G ~\~ D pG c^^j ~[~ D (^G ^vaup-^ \G ^, 

where the first identity is the usual Bianchi identity while the second one is a consequence 
of the relation g^u^o^p-iS = 9fia£iyf^-,s + Qfif^^au-yS + 9fi-y£a(3uS + QfiS^ap-tv The operator identity 
eq.(P) is new, and is one of the main results of this work. As a result of eq.@, we can 
conveniently choose a set of independent operators as eqs.(^)-(|^). For the n-th moment, 
these independent operators will mix with each other, and with the n gauge-invariant 
operators bilinear in the quark fields discussed above, under the renormalization. 

After the determination of an independent operator's basis, the renormalization is 
in principle straightforward: We follow the standard method to renormalize the local 
composite operators^]. We multiply the operators discussed above by a light-like vector 



A^. to symmetrize the Lorentz indices and to eliminate the trace terms: 

We then embed the operators Oj = A ■ R'^i, A ■ T^i, into the three-point function as, 

{0\TOjiO)A^{x)i;{ymz)\0), {0\TOjiO)A,{x)A,{y)A,{z)\0), 

and compute the one- loop corrections. We employ the Feynman gauge (a = 1) and use 
the dimensional regularization. To perform the renormalization in a consistent manner 
without subtle infrared singularities 0, we keep the quark and gluon external lines off- 
shell; in this case the EOM operators as well as the BRST-exact operators mix through 
renormalization as nonzero operators. 

One serious problem in the calculation is the mixing of many gauge-variant as well as 
BRST- variant EOM operators. As explained in Refs.[5|, P^, the gauge- variant quark EOM 
operators are given by replacing some of the covariant derivatives by the ordinary 
derivatives d'^^ in the corresponding gauge-invariant operator i?^^ In the present 

case, the BRST- variant EOM operators obtained similarly from the gluon EOM operator 
eq.(0) also participate in the mixing. However, the problem can be overcome by directly 
generalizing of the method employed in Refs. |^ Q . We introduce the vector Q!f (r = 1,2,3) 
satisfying A^Qt^ = for each external gluon line, and contract the Green's functions as 
Qt{0\T^OjA^iji^\0), Q'iQ^Q^3{0\TOjA^A^Ap\0) . This brings us three merits: Firstly, the 
tree vertices of the gauge (BRST) invariant and the gauge (BRST) variant EOM operators 
coincide. Thus, essentially only one quark (gluon) EOM operator is now involved in 
the operator mixing. Secondly, the structure of the vertices for the twist-3 operators 
is simplified extremely and the computation becomes more tractable (see eqs.([T^)-(P^ 
below). Thirdly, the three-gluon vertex of the BRST-exact operator eq.(^ vanishes for 
^2 = ^3 and Qi-Q2 = 0. Thus, we can exclude the BRST-exact operator from the gluon 
three-point functions and compute its mixing separately. 

5 



We now apply our framework to the lowest {n = 3) moment. In this case, there 
exists no three-gluon operator A • T^^. From the one-loop calculation of the one-particle- 
irreducible (IPI) three-point function using the projection by (see Fig.l), it turns out 
that the renormalization takes the following form: 
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3,E J 



where the operators with (without) the suffix "bare" are the bare (renormalized) quanti- 
ties. The renormalization constant matrix is triangular; this is consistent with the van- 
ishing physical matrix elements of the EOM operators and of the BRST-exact operators 
(see below). 




Figure 1: IPI diagrams for three-point functions (full lines: quarks; wavy lines: gluons) 



In the MS scheme, we express the renormalization constants Zij as, 

^2 



Zij — Sij -\- 
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-Xij {i, j = 1, m, B, Ep, Eg) 



' 87r2(4-L>) 

with D the space-time dimension. The relevant components of read: 

12 2 2 



X 



lEp 



X 



lEa 



-l^r, x„ 



-QCp, 



(10) 



where Cp = Cq = Nc for Nc color, and Nf is the number of quark flavors. The 

physically relevant component is only Xn. However it is inevitable to identify other 
components to get a correct result for Xu. 

Now let us discuss the relation between our approach and other ones based on the 
on-shell calculation. For this purpose, we note that eq.(|^) implies the following results 
for the IPI three-point functions including one-loop corrections (we consider the massless 
quark case): 

riJh,h,h) = Z,n{Tl^)2 + ZiEa{TlEy^\ (13) 

where Tggg{ki,p,p') and ^gggi^i, ^2, k^) are the Fourier transform of the IPI Green's func- 
tions n^(0|TA ■ Rl^A^^'^\0)ipi and fi^fi^ng(0|TA ■ Rl^A^A^Ap\0)ipi with fci,p,p'(= 
—ki —p),k2 and k-^{= —ki — ^2) the incoming (off-shell) external momenta. The three- 
point vertices appearing in eqs.([I^) and (pTSj) are given by, 

Rlsf = ^(p'r9n!,i^-p9n''ii^)t', (15) 

' ' ggg S 

TlB)^2g = -9r''e'''^''^x^2^^3u{ki-ni)k+ {cyclic perm's), (16) 

T^,Ey^l = -gr''e'''^''Ax{n2,n3Ak,■n^)h + k,,nu^2■n,)Ck3-k)) 

+ (cychc perm's) - (T-I^^'^ , (17) 

where p = A-p, etc., and "(cyclic perm's)" denotes the terms due to the cyclic permutation 
of the labels for the external lines {ki,Qi,a),{k2,^l2,b), and (^3,^3,0). We now take 
the on-shell limit of the above expressions. We can identify the vectors to be the 
polarization vectors ti{ki) of the corresponding gluons. Therefore the on-shell limit is 
realized by taking p =f = kf = ki ■ Qi = for eq.(12), and kf = k2 = k^ = ki ■ Qi = 
k2 ■ ^2 = ks ■ = for eq.(|l3D. We see {^3,3)^^^ = from eq.(|l|) and thus the 



BRST-exact operator does not contribute in this limit. On the other hand, we find 



that (i?? 



(3) 



(3) 



(3) 



gqq ^^^'^/gqq' V^^^^/ggg ®' ^ vamsfi; the EOM 

operators do contribute to the three-point functions eqs. ([T^) and ( |T3|) even in the on-shell 
hmit. 

There exist, however, the one-particle-reducible (IPR) contributions of the same order 
as the IPI ones (see Fig. 2 ). 




Figure 2: IPR diagrams for three-point functions (full lines: quarks; wavy lines: gluons) 



We denote the IPR contributions corresponding to eqs.(|T^) and (|T^) as t'^^^{ki,p,p') and 
Tggglki, ^s)" Here it should be noted that the on-shell limit of three massless particles 
corresponds to a very singular (collinear) configuration in the momentum space. Actually 
the on-shell limit of the IPR diagram (Fig. 2) becomes indefinite due to the collinear 
singularity. This technical difficulty can be avoided by calculating the "non-forward" 
matrix element of the composite operator ||^. Namely we calculate the diagram with the 
composite operators at nonzero momentum {—q = ki + p + p' = ki + k2 + k-^ ) . Then 
we "define" Tg^g and ^ggg to be the g — > limit of the matrix element after putting the 
external fields on the mass shell. 

The contributions from Fig. 2 are given by, 

rg,,{k,,p,p') = -g^,f^^U''Jp,q) + W^^{-p'-q,q)j^g^^t^ 

+ ^i,^T(ki,q)j^^,9lX, (18) 
f''gggik„k2,ks) = 1 ^ H ( ^ ^ , q) ^ ^ ^ " ^ ^ ^ ^ ^ / ^ ^ ( _ ^ ^ _ q^ ^ 3 ) ^ ^ (IQ) 
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+ (cyclic perm's), 

where V^i^iawlPi 5^2,^3) = {pi -P2)M3fl'MiM2 + (cyclic perm's) is the usual three-gluon ver- 
tex. Il'^g{p,q) and n^^^(A;, g)5''^ correspond to (0|TA ■ Rl^{q)tp{-p - q)ij{p)\0)ipi and 
(0|TA • Rl^{q)A''^'{k)A^''{-k - g)|0)ipi with the off-shell external momenta. Their ex- 
plicit expressions (divergent term) read, 

KJp,q) = ^i(,)^((p + g)(^+ ^757"^+P757"^^), (20) 



nr(fc,g) = Z,^,^{-^)e^'"^A, 



kc,9'p{{k + q)'A''-{k + qW + qn} 



(21) 



where 



Zl{q) — ZiEp , Z2(q) — ZiEq- 

Substituting eqs. (pOD -(^ID into eqs.(|l8|) and (|l^), we can go over to the on-shell limit 
with q being kept non-zero. After this procedure, we take the g ^ limit leading to 

^mi = -^^E, {RIe)^2<i " ^'""^ (^^"i?)^'^^ and f-^^ = -Z^Eo (^3>)2|^- Therefore, the 
nonzero terms coming from the EOM operators exactly cancel out in the sum of the IPI 
and the IPR contributions for the on-shell external momenta: 

^m + '^m = Zii{Rlif^\ r^OT + f^3g = 0. (22) 

The results eg. (p2D correspond to the on-shell calculations of the evolution kernel in the 
literature [0], where neither the EOM operators nor the alien operators appear. The non- 
zero contributions from the EOM operators in the on-shell limit and their cancellation 
in the sums eq.(p2D are novel features in the three-point functions, whose consideration 
is indispensable for the twist-3 operators. This phenomenon is consistent with the well- 
known theorem that the physical matrix elements of the EOM operators vanish[^], because 
both IPI and IPR diagrams contribute to the physical matrix elements. We emphasize 



that, in our demonstration, a calculation of the "physical" matrix element of the composite 
operators, in general, stays with some subtleties and dangers due to infrared singularities. 
One of the consistent method to avoid this problem from the beginning might be that all 
diagrams (IPI as well as IPR) are calculated with the composite operators at nonzero 
momentum[|]. In this case, however, we must consider additional composite operators 
which differ from the previous ones by the total derivative and this fact brings in another 
complexity^ Therefore the off-shell Green's functions are much more tractable and easy 
to calculate especially for the general moment n. 

Finally, we mention the Q^-evolution obtained from eqs.(PD-(pl]): If we neglect the 
quark mass operator A ■ .Rg^, only one operator A ■ R^ ^ contributes to the Q^-evolution 
of the physical matrix elements. The relevant renormalization constant is Zn, which gives 

I dxx'gf-\x,Q') = [^^j I dxx'gr-%x,f,'). (23) 

Here b = yA^'c — |^/, and a<j(Q^) is the running coupling constant, g^^^{x,fi'^) is the 
twist-3 contribution of g2{x,fi'^) after subtracting out the Wandzura-Wilczek piece0, 
and its moment is equal to the nucleon matrix element of A ■ i?3^(yU^) up to irrelevant 
factor. The result eq.(P^) coincides with that of Refs.0, |^, though our approach is quite 
different from those works; this fact confirms the theoretical prediction eg . (|23|) , and also 
demonstrates the efficiency of our method. We note that the term ^Nf would be absent 
from the exponent of eq.(p3l) if we consider the non-singlet casef^, |^. For Nc = 3 and 
Nf = 4, the exponent of eq.(^) is while the corresponding exponent for the non- 
singlet case is Thus the singlet channel obeys rather stronger Q^-evolution compared 
to the non-singlet one. 

In the present study, we have developed a manifestly covariant approach to investigate 
the flavor singlet part of g2{x,Q'^). We derived the operator identities which relate the 

two-particle operators with the three-particle ones. We have chosen the three-particle 
^The operator ^ is unique even if g ^ 0. 
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operators as an independent operator's basis. To identify the renormalization constants 
correctly, the off-shell Green's functions are considered. In this case, it is indispensable to 
consider the mixing of the gauge invariant operators with the EOM as well as the BRST- 
exact operators. As an application, we investigated in detail the renormalization for the 
lowest {n = 3) moment. We made contact with other methods^], and also confirmed the 
prediction of the Q^-evolution obtained by them. The off-shell Green's functions are free 
from the infrared singularity coming from the coUinear configuration as well as any other 
unphysical singularities. We believe that calculating the off-shell Green's functions is the 
safest and the most straightforward method to obtain the anomalous dimensions of the 
higher-twist operators. 

We would like to thank T.Uematsu for valuable discussions and suggestions. Y. Y 
would like to thank L. Trueman and S. Ohta for useful discussions. 
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